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Abstract 

This paper is an application of the ideas of the Born-Oppenheimer (or 
slow/fast) approximation in molecular physics and of the Isaacson (or short- 
wave) approximation in classical gravity to the canonical quantization of a 
perturbed minisuperspace model of the kind examined by Halliwell and Hawk- 
ing. Its aim is the clarification of the role of the semiclassical approximation 
and the backreaction in such a model. Approximate solutions of the quantum 
model are constructed which are not semiclassical, and semiclassical solutions 

in which the quantum perturbations are highly excited. 
04.60.+n 
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I. INTRODUCTION 



A. Semiclasssical gravity and perturbed minisuperspace 



There is a well-explored and consistent theory of quantized (free) fields propagating on a 
fixed classical spacetime [|lj . In particular, one may construct the expectation value of their 
stress-energy tensor. Although it is formally infinite, there is a variety of regularization 
procedures which give the same and often sensible result. A simple example of such a result 
has been experimentally verified in the form of the Casimir effect. A more indirect link of 
quantum field theory in curved spacetime to physical reality is the possibility of explaining 
early universe density fluctuations as quantum fluctuations. 

It is tempting to insert the quantum stress-energy tensor into the classical Einstein 
equations which govern the background spacetime. In the Heisenberg picture we could 
formally write 



where g is the classical metric, q the quantum field operators and \if) > their quantum state. 
To make the theory a consistent approximation, one must include among the quantum fields 
the linearized perturbations of the metric itself. Higher order effects from the quantum na- 
ture of the gravitational field should in contrast be suppressed by a factor of the background 
curvature scale over the Planck scale. It has been argued that this scheme, called semiclas- 
sical gravity, is a better approximation to the presumed quantum theory of gravity than 
classical general relativity if the curvature scale of the background spacetime is much larger 
than the Planck length. The question if it really is must be considered open in the absence 
of both experiments and a quantum theory of gravity. The effective equations of motion for 
the background metric, taking into account the backreaction of the quantum fields, contain 
higher than second time derivatives, which give rise to qualitatively different new solutions. 
It has been argued that these are artefacts of the approximation to full quantum gravity 
and should be excluded in a consistent way f2[. 




(1) 
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A tentative model for quantum gravity in which the derivation of semiclassical gravity as 
an approximation to a more general theory can be attempted is the canonical quantization of 
a minisuperspace model with generic but linearized perturbations. One such model has been 
given by Halliwell and Hawking ||, where the minisuperspace model is a closed Friedmann 
universe driven by a self-interacting real scalar field which is homogeneous on the surfaces 
of homogeneity. 

The Halliwell-Hawking model does not seem to be totally consistent: It has been pointed 
out [f| that Halliwell and Hawking miss some components of the Einstein equation to a self- 
consistent order by restricting their ansatz for the background metric before variation. A 
related problem is the presence of linearization instabilities due to the symmetries of the 
spatial hypersurfaces |]5]|| . In this paper we shall start from the quantized Halliwell-Hawking 
model as it stands in spite of these problems. Our reason for this neglect is simply that 
their resolution will require the imposition of a finite number of additional constraints on 
the wave function, in addition to the infinity that is already there. Although essential for 
consistency, these should not change the physical content of the theory, in which we are 
interested here. 

Our aim is to reconsider the role of the semiclassical backreaction by a method which 
applies to the Halliwell-Hawking model, but which, in its physical meaning, should generalize 
to other perturbed minisuperspace models. 

Halliwell J?|,[8| and other authors PHTT] have given a derivation of the equations of semi- 



classical gravity which is based on making a single-factor WKB ansatz for the wave function 
as in reference 0. As we shall see directly, the backreaction term found this way is small 
by virtue of the WKB approximation itself, and the theory obtained therefore less general 
than semiclassical gravity. Semiclassical gravity by naive insertion of the quantum stress 
tensor into the Einstein equations would in contrast also allow for highly excited states of 
the quantum fields which contribute a large part or all of the total source term of the Ein- 
stein equations. It is this possibility which we examine here, based on the Wheeler-deWitt 
equation of reference |§. For simplicity we consider only pure gravity, so that the Fried- 



mann universe will be driven by gravitational waves only. Classically, this possibility has 
been examined a long time ago ||12|| . The approximation method we shall use is related both 
to the Born-Oppenheimer approximation in molecular physics [[HJ and to the short-wave 
approximation of Isaacson for gravitational waves ||I4[. In the remainder of this introduc- 
tion we review the previous approach to the backreaction derivation and then the relevant 
features of the Born-Oppenheimer and Isaacson approximations. 



B. The Halliwell-Hawking approach to the backreaction 

Let us consider a model in canonical gravity where the Hamiltonian constraint is of the 
form 

H(Q,P, q, p) = H (Q, P) + H 2 (Q, P,q,p) = 0, (2) 

where H 2 is homogeneous of order 2 in the p and q. In reviewing the general way in which 
Halliwell and Hawking find an approximate solution to the corresponding Wheeler-deWitt 
equation, we can limit ourselves to a toy model which will display the features relevant 
for our argument. For simplicity of notation alone we assume therefore that the quantum 
operator version of Hq is simply 



(We have assumed a positive definite kinetic energy term. In gravity it is indefinite.) With 
the ansatz ^(Q, q) = e tS ^ip(Q, q) we then obtain 

l -S 12 ^ - iS't// - ^" - l -S"ij + Vip + H 2 ij = 0, (4) 

where a prime denotes d/dQ or d/dQ. Dividing by ip we can separate this equation into 

l S a-*S» + V + E (Q) = 0, (5) 
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ip = 0. (6) 
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Here E(Q) is an arbitrary separation function. It has no direct physical significance, as its 
choice does not affect the product e tS ip, but an appropriate choice can help in the interpre- 
tation of the two separate equations of motion. We now assume that both e tS and the total 
wave function e lS ip are of the WKB form in the variable(s) Q, that is we assume 

\S' 2 \ > \S"\, (7) 

> |(lnV)'l- (8) 
To leading order in this approximation we obtain the Hamilton- Jacobi equation for S(Q), 

l S ' 2 + V + E(Q) = 0. (9) 

By the argument S' — P — dL/dQ — Q (a dot denotes d/dt or d/dt), we can identify 
S'd/dQ with d/dt in the WKB case. To leading order we then obtain a time- dependent 
Schroedinger equation 



-i^ + H 2 - E(Q) 



4. = 0, where | s S> ±. (10) 



One might be tempted to set E(Q) =< il)\H 2 ip > in (0) in order to obtain the Hamilton- 
Jacobi form of the semiclassical backreaction ([I]) in the Einstein equations in this way, as 
was once proposed by Hartle [ 16| . But in (|10|) this leads to < if} \ — id/dt ip >= 0, so that the 



other part of the semiclasscial gravity scheme, the time-dependent Schroedinger equation 
describing the evolution of the quantum fields, is lost. If we want to obtain the correct 
time-dependent Schroedinger equation for the gravitons, i. e. one where the Hamiltonian 
is the usual Hamiltonian of linearized gravitational waves in the given background Q(t), 
we must set E(Q) = 0. But then the background evolves without a backreaction from the 
perturbations. In order to obtain semiclassical gravity as it is commonly understood, one 
would have to set Q(Q) equal to in ([Uf), and equal to < ip\H 2 ip > in @, a question of 
having the cake and eating it. 

A subtle argument for the presence of a backreaction in spite of this dilemma has been 
given by Halliwell 0. He determines what the background spacetime really is, not from 



equation (||), but by a Wigner function analysis of the total wave function. (The convention 
for E(Q) is then irrelevant.) Following Halliwell, one may assume that the perturbations q 
are for some reason unobservable in principle. Predictions should then be made, not from 
the wave function ty(Q,q), but from the density matrix which is obtained by tracing over 
the perturbations, p(Q, Q') = J ^(Q, q)* ^(Q f , q) dq. The Wigner function derived from this 
density matrix may then be positive definite and may be peaked around a classical trajectory 
in phase space (Q,P). Halliwell showed that this trajectory is not given by P{Q) = S'(Q), 
but rather by 

P(Q)=S'(Q)+<4,\-i^>. (11) 

Putting this into the Hamilton- Jacobi equation (|), using (|1(J) and treating the second term 
in (|TT|) as small, one obtains 

^P 2 + V+ < tfj\H 2 ^j >= 0. (12) 

(As predicted, the choice of E(Q) does not matter.) As a way of finding the correct physics in 
spite of the appearances of the equations of motion (p|JT0|), this is an impressive argument. It 
is limited in two ways. By the assumption \S'\ 3> |(ln^)'| which went into its derivation, the 
backreaction term derived in this way is necessarily small and cannot qualitatively change 
the background. This limitation is also apparent in the fact that the equations one has to 
solve are still (^) and (|10D , and not those of semiclassical gravity. 

Secondly, only those quantum perturbations contribute to the backreaction which are by 
assumption unobservable. There are good tentative reasons for this assumption. In a closed 
universe which has undergone inflation, for example, one would not be able to observe the 
quantum state of the perturbations completely because of the presence of particle horizons 
H . But why should one only ever be able to observe the gravitational backreaction of those 
quantum fields that are "out of sight"? 

The limitation of the backreaction effect to a small perturbation of the background 
trajectory is even less satisfactory. One should expect to find a quantum equivalent of 
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the classical situation where the self-gravity of gravitational waves substantially curves the 
background spacetime on which they move. Examples of such spacetimes which have been 
examined are a nearly static spherically symmetric concentration of gravitational waves (a 
gravitational geon [fT7||), or a Friedmann universe closed by and driven by gravitational waves 

0- 



C. The classical Isaacson approximation 



A key ingredient for both these cases was supplied by Isaacson [14j (see also @]). He 



considered a one-parameter family of solutions of the vacuum Einstein equations of the form 

9ab{x; e) = 'jab(x) + h ab (x; e) (13) 

and assumed that the family of perturbations h ab (x;e) is such that their amplitude scales 
like e, but their "wavelength" like e -1 . Formally, one assumes that there is a coordinate 
system in which 

h a b(x; e) = 0(e), h abfi (x] e) = 0(1), h ab , cd (x] e) = 0(e _1 ), (14) 

where a comma denotes a partial (coordinate) derivative. ^ ab and its derivatives are consid- 
ered as 0(1). Let the Ricci tensor of g ab be expanded in powers of h ab : 

Rab(9ab, gab,c, g a b,cd) = R%(l) + Kbil ') h ) + R lb(l h ) + ■ • • = °> ( 15 ) 

where R\ b is linear in /i a fc, R 2 ab quadratic, and so on. Here R\ b is a contraction (with r y ab ) of 
V a V b h cd , where V ' a lbc = 0, and R 2 ab is a contraction of V a V b h c dh e f . Inserting the ansatz 
(|i~3|) into the vacuum Einstein equations, and separating powers of e, one obtains to leading 
and next order 

R l ab(l) = (16) 
R°ab(l) = -< Rlbil, h) >average (17) 
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(Strictly speaking, the connection terms in R\ b and R 2 ab are of a lower order than the partial 
derivative terms and do not appear in these equations). The brackets in ( |TTD (and here 
only) refer to a suitable spacetime averaging. The right-hand side of ( |T7| ) is an effective 
stress tensor of the perturbations which drives the background. Equation ([RJ) is a linear 
wave equation for the perturbations propagating on the background. These equations have 
to be solved together in a self-consistent way, and as such are not part of a perturbation 
scheme in which each order is solved before and independently of the higher ones. The 
difference from conventional perturbation theory is also seen from the fact that the limit 
e — > is singular. 



D. The Born-Oppenheimer approximation 

It is worth reviewing the Born-Opppenheimer approximation in a notation which is 
sufficiently abstract to be suggestive both of its original context of molecule physics and 
of the context it will be used in in this paper. Let us assume that we are dealing with a 
Hamiltonian of the form 

H(Q,P,q,p) = H (Q,P)+H 2 (Q,q,p). (18) 

The single but crucial difference to (0) is that H 2 depends only on Q, but not on the conjugate 
momentum (or momenta) P. We want to solve the time-independent Schroedinger equation 

H(Q, -i-^, q , -i-^MQ, q) = e*{Q, q), (is) 

where E is a given constant. (The case of canonical quantum gravity differs from this only 
by the presence of additional constraints on the wave function, and by the restriction E — 0.) 
We use the fact that H 2 does not contain any Q-derivatives to pose the eigenvalue problem 

H 2 (Q, q, -^)MQ, q) = E v {Q)f v {Q, q) (20) 

for the complete sets f v and E u , for each value of the parameter Q. We choose the /„ to be 
orthonormalized under the scalar product 
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</„!/,>= / rj.dq. (21) 
(They are of course already orthogonal for E v ^ E^.) We now make the ansatz 

*(Q,?) = EWW.?) ( 22 ) 

V 

and solve for the Let us again assume, for simplicity of notation, that Ho is of the form 
(0). (It would be straightforward to write down what follows for H being an arbitrary second 
order derivative operator in more than one variable. This general case includes canonical 
gravity.) We can simplify the result obtained from substituting (|22|) into ( [L9|) by splitting 
it into components with respect to the basis f u using the scalar product (|21[) . We obtain an 
infinite number (labelled by v) of coupled equations in the minisuperspace variable(s) Q: 

(H + E V {Q) - E) *„(Q) = E < Ul-^U > < U\^qU > (23) 

In the molecular physics application, where Q are the nucleus positions and q the electron 
positions, it can be shown that the right-hand side is small by a factor of the square root 
of the ratio of the electron mass over the proton mass. The argument itself does not seem 
capable of extension to gravity, although it has been argued that the Planck mass must be 
involved in finding a small number in the problem. This is obviously inapplicable to the 
case of pure gravity, where the Planck mass is the only fundamental scale. In our particular 
example of a perturbed Friedmann universe the suppression factor will turn out to be the 
ratio of the length scales of the metric perturbations over the Hubble scale of the background 
metric, a very intuitive result given the usual language in which the q are called the fast 
variables and the Q the slow variables. Neglecting the right-hand side coupling terms in 
( p3|) , and assuming E = 0, we can suggestively write it as 

(H +<f u \H 2 f u >)V u (Q) = (24) 

Here we seem to have obtained the backreaction equation we were looking for in a direct 
manner and without making any arbitrary choices. However, we are not really solving two 
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equations self-consistently, in the manner of (|16|) and flT7|). Instead we have to do more, 
namely solve ( p0|) for all values of Q. There is no classical equivalent of this. 

In section II we describe the Hamiltonian formulation of a perturbed minisuperspace 
model - that of Halliwell and Hawking without matter - which is is completely described by 
a Hamiltonian constraint of the form (pj|). We make the first steps towards the Born- 
Oppenheimer approximation by solving the fast-part eigenvalue problem (|20| ) and then 
writing down the slow-part equation (p3|) . In section III we find the Born-Oppenheimer 
approximation to that model, by finding the regime in which the right-hand side of (p3|) can 
be neglected. This turns out to be the Isaacson approximation. In section IV we recover 
semiclassical gravity by adding the WKB assumption. We compare our results with previ- 
ous work. Section V resumes our results and raises some questions on the validity of our 
starting point (equation ( [26]) below) and the posssibility of generalizing our treatment to 
other perturbed minisuperspace models. 

II. BORN-OPPENHEIMER TREATMENT OF A PERTURBED 
MINISUPERSPACE MODEL 

Perhaps the simplest physically complete perturbative minisuperspace model is that of 
a Friedmann universe without matter, but with generic (small) perturbations of the metric. 
For a closed (K = 1) Friedmann universe this is a special case of the model of Halliwell and 
Hawking [Q. There is a single Hamiltonian constraint of the form (@) - where if 2 is the 
integral over all space of the quadratic Hamiltonian constraint, and Ho is the minisuperspace 
Hamiltonian constraint - and an infinity of constraints (the linearized Hamiltonian and 
momentum constraints at each space point) which are homogenous of order 1 in the p and 
q. (There should also be a finite number of additional constraints which are homogeneous 
of order 2, but as stated in the introduction, with Halliwell and Hawking we neglect these 
here.) 

A canonical transformation that both solves the linearized Hamiltonian and momentum 
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constraints and brings the remaining part of the Hamiltonian constraint into the "Born- 
Oppenheimer" form (|IBD has been given by Wada JOJ. The wavefunction then depends 
only on some of the new variables, and is subject only to a single Hamiltonian constraint. 
The remaining degrees of freedom in this particular model are the amplitudes of transverse 
traceless perturbations of the 3-metric, d n , where n labels the space dependence (the nth 
tensor harmonic on the 3-sphere), and the minisuperspace variable a, which is equal to a, 
the logarithm of the scale factor, up to terms quadratic in the perturbations. 

Although physically interesting in its own right, this model will here also serve as a toy 
model for a more general perturbative minisuperspace ansatz. Therefore we use the previous 
notation Q for the minisuperspace variables and q = {q 3 , g 4 , . . . q n , . . . } for the perturbation 
variables, (n is an infinite discrete label which arises from splitting the perturbations into 
"Fourier" components with respect to their spatial dependence.) In the following Q is 
therefore the same as a and q n is the same as d n of reference ||19|| , where their precise 



definitions can be found. The definition of the perturbed spacetime metric in terms of these 
variables is given in [[J . For understanding the physics of this particular model it suffices to 
know that the background metric is (now using our notation) 

ds 2 = -N 2 dt 2 + e 2Q dQ 2 3 , (25) 

where dVl\ is the round metric on the unit three-sphere, and where N and Q are arbitrary 
functions of t, and that q n is the amplitude of a transverse traceless perturbation of the 
3-metric of comoving wavelength ~ 1/n. The total Hamiltonian of the perturbed system, 
in the remaining degrees of freedom given by Wada (but in our notation) is 



H = Nne 



-3Q 



1 1 00 

l -P 2 -e^ + W(pl + ^ 2 -l)e^q 2 n 



2 2 

Z Z n=3 



(26) 



It is constrained to vanish. We shall find a family of approximate solutions - in the Isaacson 
limit - of the corresponding Wheeler-deWitt equation 

'1 d 2 An 1 g ( d 2 

n=3 



2 9Q* h^ 2 - 1 )^ 



*(Q,g) = 0. (27) 
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We have made an ad-hoc choice of factor ordering. 

As a first step, as yet without an approximation, we split the Wheeler-de-Witt equation 
in an infinity of variables into an infinity of equations in only the minisuperspace variables, 
here only Q. This is precisely the ansatz that prepares the Born-Oppenheimer approximation 
in molecule physics. We make the ansatz 

oo 

9(Q, q) = E MQ) II U ((™ 2 - l) 1/4 e Q gn) . (28) 

v n=3 

Here v = {z/ 3 , i/ 4 , z/ 5 , v n , ...} is an infinite tupel of numbers which have range v = 0,1, 2, .... 
f v {x) is, for each value of v, a given function of its argument x, namely the normalized 
energy eigenfunction of the unit frequency harmonic oscillator, defined by (1/2) (— d 2 /dx 2 + 
x 2 )f u (x) = iy + l/2)f u (x) and J f 2 dx = 1. In this definition we have already incorporated 
the solution of the eigenvalue equation (|20|) for H 2 given by fl27j). As the f v form a basis 



of square-integrable functions on the real line, our ansatz (|28| ) is generic. We note that 
d/dQf v = xd/dxf v where x is the formal argument of f Vn , namely (n 2 — l) l / A e Q q n . We 
introduce the operators a = (x + d/dx)/y2, a) — (x — d/dx)/y/2 and N = a^a, with the 
well-known properties af v = \fv f v -i, & fv = \/v + \ fv+i, afo = 0, and Nf u = vf v . Now 
we follow the Born-Oppenheimer procedure of separating the <Ev by the orthonormality of 
the f v . For this purpose we use the identity xd/dx = (l/2)(a 2 — a^ 2 — 1) and its square 
(xd/dx) 2 = (l/4)(a 4 + a t4 - 2a 2 + 2at 2 - 2iV 2 — 2N — 1). The final result is 

(1 d 2 „„ 1. 



- - e 4 Q + e 2 « Y>n + ) *p = ~\ £[^>„ +2 - 



\2dQ 2 " V 2 / 2 

~o £|$*K-H + ®v\v n -4 - 2$ff\ Vn+2 + 2$j^_ 2 - (2z/ 2 + U n + l)$ff 

8 n 

[®v\v n +2,v m +2 + $ff\v n -2,v m -2 - 2®v\v n -2,v m +2 ~ 2®v\v n +2 + 2 < l;7|fc Vl _2 + $ 



(29) 



Again, a prime denotes d/dQ or d/dQ. To keep this expression readable, we have introduced 
a shorthand notation that implies certain factors arising from the Bose-Einstein statistics of 
the excitations of the perturbation modes q n . It is most easily explained by example: 
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u\u n -4 z 



v \Jn v - \\fn v - 2\pn v - 3$ 



{ ^3,^4, ■■■,V„-l,V n ~4,V n +l,...}l 



d 



(30) 
(31) 



The left-hand side of equation ( |29| ) contains an infinity of the form J2™=3 y/n 2 — 1, and 
the right-hand side several infinities of the form J2^=3 1- We must eliminate these if the 
equation is to have more than formal meaning. A slightly more consistent procedure than 
to just cross out these terms is to normal-order all occurences of a and a\ i. e. to put a to 
the right of a) in all products. In H 2 we obtain N instead of N + 1/2. This is in analogy 
to standard regularization procedure for free quantum field theory in flat space, xd/dx 
becomes : xd/dx := (1/2) (a 2 — a^ 2 ). Replacing (xd/dx) 2 by : (xd/dx) 2 : does not get rid of 
all infinities. Instead we choose 

: (: xd/dx :) 2 := (l/4)(a 4 + a t4 - 2N(N - 1)). (32) 

In support of this choice one might argue that xd/dx should always be normal-ordered as a 
block as in (0), because it arises from the expression df v /dQ. With that choice of normal 
ordering we obtain 

/ 1 ^ e*Q + e 2 « £ u n y/tf=l) d>, = -~ £ [% K+2 - % 



2dQ 2 Z^"nVn ij^ v — ^v\u n -2 

8 n L 

_ o X! X! [^K+2,i/ m +2 + $^-2,1^-2 - 2<IV| i , n _2,i/ m +2 



(33) 



In the following we use this equation as our starting point. We shall refer to the v n as "oc- 
cupation numbers" by analogy with the states of a set of harmonic oscillators. Furthermore 



we shall by the same analogy refer to the first line of the right-hand side of (p3|) as the 
"two-particle terms" and to the remainder as the "four-particle terms" . 
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III. QUANTUM VERSION OF THE ISAACSON EXPANSION IS THE 
BORN-OPPENHEIMER APPROXIMATION 



We now introduce the Isaacson expansion as a means of making the right-hand side 
of (|33|) small compared to the last term on the left-hand side, thus completing a Born- 
Oppenheimer approximation scheme. (To post-Born-Oppenheimer order we shall recover 
semiclassical gravity.) With Isaacson we consider a family of approximate solutions labeled 
by a small parameter e, such that the large-scale metric, in our case the Friedmann back- 
ground, is independent of e. Both the amplitude and the wavelength of the perturbations 
are chosen to scale like e, so that their effective energy-momentum tensor, suitably aver- 
aged, is independent of e and can act as the source for the background curvature. In our 
model, defined by the Hamiltonian (p6|), we identify as a guidance principle the perturbation 
wavelength with 1/n and the classical amplitude with the occupation number v n . The total 
energy of the perturbations is then proportional to 



This is of course the term appearing on the left-hand side of (|3~3"D which we would like to 
retain as a backreaction term after the model of (p4|). To make this expression (approxi- 
mately) independent of e, we may consider the one-parameter family of occupation numbers 
(classically, amplitudes) 



The explicit e-dependence in the second argument of the function fi(n; e) is necessary to 
make the perturbation energy precisely e-independent with this ansatz. Isaacson however 
uses the fact that for small e this dependence becomes weak, so that for sufficiently small e 




(34) 




(35) 



the function /i(n, 0) characterizes a whole family of approximate solutions vffi by the scaling 
relation ~ e 2 /i(en, 0). 



The equivalent in the quantum theory as defined by (^) is a one-parameter family of 
wave functions $^ (Q)- We define it as 
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4 £) (Q) = Fp(Q; e) where /x n = e" 2 z/ £ -i n . (36) 

To keep both the value and the suffix of fi n integer one must restrict the range of e to 

e = N' 1 , with N any positive integer. A similar problem, that the family of approximate 

solutions cannot be continuous in e, would arise also in the classical Isaacson approach to 

a closed universe. It therefore does not invalidate the existence of the approximation. In 

the case of open spatial hypersurfaces, the suffix n becomes continuous anyway, and one 

may then simply consider F(fT) as a smooth functional on a continuous field. Alternatively 

one may set e = 1 at the end of a formal expansion in e. Then one no longer considers a 

one-parameter family of solution, but rather a single solution whose occupation numbers are 

nonvanishing for high frequencies only, and then small. 

By substituting the ansatz (^) into the equations of motion for $j7, and then changing 

variable from z7 to /2 and summation index from n to en, we obtain an equation of motion 

for F that contains explicit powers of e in its coefficients: 

(1 d 2 
\2dQ 2 

= ~\t 2 Yl[\[^n + 2e- 2 ^Jfi en + e- 2 i^| M£n+2e - 2 (Q; e) - ^JI^^J fi en - e- 2 F'^ cn _ 2e . 2 {Q; e 

Z n 

-l^TX^^n + 4e" 2 v //i en + 3e- 2 ^/fi en + 2e~ 2 ^ en + e^F^+^Q; e)] - 

° n ° n m^n 



- e 4 Q + e 2Q £ u n y/tf=l) F n (Q- e) 

n J 



We have not written out all terms. The remainder can easily be reconstructed from (j33"[). 
The important point is that all "two-particle" terms of the right-hand side are multiplied 
by a factor of e 2 and all "four-particle terms" by a factor of e 4 . These factors arise from the 
Bose statistics prefactors, which we have written out here for this reason. 

We see from the presence of these explicit factors of e in ( p7| ) that it was necessary to 
have given F an explicit dependence on e (as well as the dependence of its formal argument 
ft on e). But by analogy with the work of Isaacson, we expect that solutions for different 
values of e are related by a simple scaling of their arguments, as e becomes small. Expecting 
the same physical phenomenon in the quantum theory, we have already incorporated that 
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scaling behaviour into (|36|). We now express our expectation that the explicit e- dependence 
of F disappears as e — > by expanding it as 

F${Q\ e) = X p{Q) P^ 2 Q)a^Q)... (38) 

The final form of the equations of motion with our ansatz is then, after separating powers 
of e and then setting e = 1, 

(^2 + e2Q £ "nV^T - e 4 «) Xi7 (Q) = 0, (39) 

Xu(Q) ^gM<5) = £ (x^ n+2 ( ( 5) - X*K- 2 (<3) P*K-2(Q)) , etc. (40) 

The leading order equation ( |39"D is of course just (|37D with e set equal to zero. As such 
it is the Born-Oppenheimer approximation to the quantum equations of motion. From the 
analogy with the classical Isaacson expansion we expect that it describes the motion of the 
background driven by the averaged perturbations. The next order (f40"D should describe "the 
motion of the perturbations on that background" . In particular, if the solution to ([!9|) is 
of the WKB form, thus implying a classical background Q(t), we should expect (pUf ) to 
describe a free quantum field theory on the curved classical spacetime given by (p5|). The 
following order, a first-order linear differential equation for a with coefficients depending on 
X and p, should describe higher order corrections that we do not expect to understand in 
this intuitive sense, and we have therefore not written it out. It may be worth clarifying 
that, although the physical significance of e here is the same as in ([14]) , mathematically it 
leads to a perturbation expansion as in (||[10]) rather than a self-consistent field expansion 
as in (0,0). 

Before we examine the semiclasscical gravity interpretation of (|39|), fl40"D, we note here 
that this set of equations has solutions which are definitely not semiclassical. We have 
therefore found approximate solutions of the pure gravity Halliwell-Hawking model which 
are not semiclassical, as promised in the abstract. In other words, it does not actually make 
a difference to solving the perturbation equations if the background is behaving classically 
or not, because we need not look for a time-dependent Schroedinger equation. 
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IV. WKB APPROXIMATION AND SEMICLASSICAL GRAVITY 



We now compare our results with the corresponding quantum field theory of perturba- 
tions on a fixed background spacetime, with a given Q(t) (for a given choice of No(t)) in 
the background metric fl25|). As a first step we must formulate precisely what we mean by 
quantum field theory in a curved spacetime in the context of our model of a Friedmann 
universe with gravitational wave perturbations. The Friedmann universe (^) is to be the 
classical background, and the metric perturbations are to be the quantum fields. As we 
want to make contact with quantum cosmology, we choose the Schroedinger picture. 

Instead of a Wheeler-deWitt equation the perturbations obey the time-dependent 
Schroedinger equation 



4 + >o(t)e" 3 ^ ± (-i- 2 + (n 2 - 1) £ 



tt(*,g) = 0. (41) 



dt 2 - ' f 3 V ^ 2 V 

We again split the wave function into harmonic oscillator components 

oo 

q) = £ Mt) II ((n 2 - l) 1/4 e QW g.) (42) 

v n=3 

with the same notational conventions as before. We may then apply the same procedure 
of identifying the action of d/dQ on the f Vn with the action of annihilation and creation 
operators. The only difference is that d/dQ arises here from —id/dt acting on the Q(t) 
inside the f Un . Therefore only first Q-derivatives appear, and the equations of motion for the 
coefficients are much simpler. We use normal ordering again as a regularization procedure. 
As there are no terms quartic in a and a', it is unambiguous. We obtain 

In these equations we recognize the well-known fact that "particles" are created in pairs 
(due to the symmetry of the background spacetime) and mainly in spacetime regions where 
the curvature scale is less than or comparable to their Compton wavelength. This is not 
the place to review definitions of the particle concept in curved spacetime. It will suffice to 
consider a simple situation in which Q(t) is constant, then changes, then is constant again. 
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In the regions of constant Q we have a well-defined notion of ground state and particles. 
The ground state is that in which all vanish apart from $g. If we start with this state at 
an early time, then after the period of nonvanishing Q, at late times, the other components 

will no longer vanish. The components with v n = 2,4, .. . for a given n will be excited 
by being strongly coupled to the z7 = component, if \dQ/dt\ > N e~ Q n during that period. 
But e^/n is the physical wavelength of the perturbation labelled by n, and Nq 1 dQ/dt the 
Hubble constant (with respect to proper time). Therefore the criterium that particles with 
comoving wavelength 1/n be produced abundantly is just that their physical wavelength be 
greater than the Hubble distance. 

It is illuminating to separate off the trivial part of the time evolution of the (that is 
the part which is independent of their initial values): 

= pff(t) exp (-i VnV^l f N (t')e- Q W dA (44) 



(pff defined here is not the same as in equation fl38|), but we have chosen this notation 
consciously, as the two definitions will be seen to be related.) The resulting equation for the 
pa is now real: 

j t Pu{t) = —-Q E (A?K+2 - &K-a) (45) 

One way of looking at this equation is as an interaction picture of quantum field theory: The 
trivial part of the dynamics has been suppressed and only the interaction part is described 
explicitly. In flat space quantum field theory this is the nonlinear part of the equation of 
motion, but in curved space also the interaction with the background curvature, even for 
linear fields. 



Although equation (j"5|) is real, it still describes the entire evolution of the quantum 



fields q on the background spacetime ( pip - the multiplication by the factor ( f44|) is merely 
algebraic. But we have now formulated the time-dependent Schroedinger equation in a 
way which evokes Hartle's original attempt at finding the backreaction in an unexpected 
manner. In particular, if we define ip(q,t) = Y^vPvfv-, we fi n d < ip\d/dtip >= 0. The energy 



of the perturbation has been hidden away in equation (fP|), and we can always reconstruct 
it trivially. 

In order to establish the relation of (Jfijj) to the full Wheeler-deWitt equation, we can 
use Q(t) locally as the independent variable instead of t, by replacing d/dt by (dQ/dt)d/dQ. 
The equation for we obtain is then very similar to equation 

^qPv(Q) = ~\ - £*K-a) (46) 

This resemblance is rather surprising, as we have not yet introduced a semiclassical 
approximation, only the Isaacson expansion. What is missing to recover (f43l) from (|40|)? 
The first requirement is that x'u — Xp\ Vn ±2 f° r a ^ n ' so that these factors cancel in fldOl) , 
allowing us to go from there to (46). Clearly this is possible for a certain interval of Q given 



appropriate initial conditions, because the equations (|39|) for and Xu\u n ±2 have nearly the 
same coefficients: The difference is only the energy of two gravitons compared to the energy 
of all gravitons in the universe. Our first requirement is therefore that Xu{Q) be a smooth 
function not only of Q, but also of V. 

The second requirement , for going from fl4"6|) to ( [45] ) , is not surprising: Clearly dQ / dt is 
not defined for arbitrary Xu{Q)t but only by the Banks || 1 5|] - Halliwell- Hawking [[5] approx- 



imation. To obtain it, we must assume that the xv are of WKB form. Finally, we obtain 
(flU) from fl4"5|) by the trivial step (§4). 

What is the equation of motion for the semiclassical background? We substitute the 
WKB ansatz Xv — e~ lS ^ into the minisuperspace equation (|39D and in the WKB approxi- 
mation \S'\ 2 ^> \S"\ approximate x" by —S' 2 x and interpret S' as the classical momentum 
P. From the Hamiltonian (p6|) we read off Q = — N e~ 3 ®P and obtain finally the classical 



energy constraint 

2 



+e- 2 « = e-^E^v^T. (47) 

As e®® is the Friedmann scale factor, and Q/Nq the Hubble constant, this is just the 
Friedmann equation for a closed universe in the presence of a radiation fluid, whose density 
scales as (scale factor) _4 . 
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In rederiving semiclassical gravity we have not quite arrived at (jl]). From a naive inter- 
pretation of that equation one might expect that the backreaction term of a mixed state is 
just the sum of the backreaction terms of its pure-state components. This seems to be true 
in our model only if the components have approximately the same energy. In other words 
we cannot have the background driven by an effective perturbation stress tensor which is an 
average of the stress tensors of two different pure quantum states, unless these stress tensors 
are nearly the same. Presumably such a situation would not make experimental sense. 

It should also be noted that we have found neither negative energy in the effective stress 
tensor of the quantum fields - it is just that of a radiation fluid - nor fourth-order effective 
background equations. This may be due to our choice of model to start from. 

V. CONCLUSIONS 

We have found a more general way of (approximately) solving the Wheeler-deWitt equa- 
tion of the Halliwell-Hawking model than the one making a WKB ansatz for the back- 
ground variables. Our approximation requires that the system splits into a slow and fast 
part, namely the Friedmann background metric and short-wave metric perturbations, with 
intermediate wavelengths in or close to their ground states. Physically, this assumption 
corresponds to the Isaacson approximation in classical gravity, mathematically to the Born- 
Oppenheimer approximation in quantum mechanics. 

Under this condition the Wheeler-deWitt equation splits into a system of uncoupled 
minisuperspace Wheeler-deWitt equations for the Friedmann background and a system of 
coupled first order equations for the perturbations, with coefficients depending on the solu- 
tions of the background equations. The essential difference to the Halliwell-Hawking solution 
scheme is that one does not have to assume that the background is in a WKB state to obtain 
this split. 

The WKB form of the background can be chosen as a second, and independent assump- 
tion, which allows to interpret the minisuperspace WdW equations as a Hamilton- Jacobi 
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equations for a classical background driven by a classical radiation fluid, and the system of 
perturbation equations as a time-dependent Schroedinger equation for QFT on that clas- 
sical background. In other words, the WKB approximation works in the same way in our 
scheme as in the Halli well- Hawking scheme, but the class of semiclassical solutions we ob- 
tain is more general: The Hamilton- Jacobi equation (the Friedmann equation in our model) 
contains the energy of a classical radiation field. This term, the semiclassical backreaction 
of the quantum metric perturbations, can be large, i. e. the quantum fields can be highly 
excited. 

The most important remaining question is whether a Born-Oppenheimer expansion can 
be applied to any other perturbed minisuperspace models. We do not know why the canon- 



ical transformation described by Wada |19| manages to solve all the linearized constraints 
of the Halliwell-Hawking model restricted to pure gravity and at the same time to bring 
the Hamiltonian into a form that allows our Born-Oppenheimer treatment. The same kind 
of transformation, applied by Shirai and Wada [5I| to the full model including scalar field 
matter, is not so successful: The linearized Hamiltonian constraint is not solved, and the 
zero-plus-second-order Hamiltonian is not of Born-Oppenheimer form in the remaining de- 
grees of freedom. There is no obvious reason for this failure. 

Wada's treatment is hard to follow in geometrical terms, as he works from the Hamil- 
tonian of Halliwell and Hawking, which is complicated for two different reasons. Firstly, 
Halliwell and Hawking split all perturbations into harmonics on S 3 from the start. Sec- 
ondly, they truncate the action and then calculate the Hamiltonian. One would expect 
more insight from truncating the full Hamiltonian directly, a procedure which may be in- 
equivalent if the background is not itself a solution. In truncating the full Hamiltonian 
directly, the source of terms of the form f{Q)Ppq and f(Q)PPqq in H 2 is clear: They 
arise from expanding the "supermetric times momentum times momentum" part of the full 
Hamiltonian constraint. Still, one may able to see for the general case of a homogeneous 
background if a canonical transformation cannot be found which gets rid of these terms. 
At least one should be able to formulate Wada's transformation in geometrical terms and 

21 



explain why it does not generalize, if indeed it does not. In such a treatment it would be 
necessary and probably illuminating to correct the inconsistencies of the Halliwell-Hawking 
model we mentioned in the introduction. 

After this work was completed, reference was brought to the author's attention. 
There a general, exact solution of the quantum problem we have considered here is given. If 
one makes the canonical transformation q n = (exp Q)q n , Q = Q + \ J2 n Ini the Hamiltonian 
(^H) transforms into a sum of simple harmonic oscillator Hamiltonians, the one representing 
exp Q with a negative sign. (The second part of this transformation is not made explicit in 
||21|| .) One can therefore give the general exact solution of the corresponding Wheeler-deWitt 
equation in terms of products of harmonic oscillator eigenf unctions. 

The simple expression for the transformed Hamiltonian arises in the approximation in 
which terms of order J2nQn are kept, but terms of order (J2nQn) 2 are neglected. This is the 
same approximation that was used in truncating the full Hamiltonian in the first place, i. 
e. to keep terms to quadratic order in the metric perturbations. It breaks down towards 
the initial or the final singularity or both of any classical solution, where the perturbations 
become large and must be treated nonlinearly. 

The difference between exp Q and exp Q is of order J2 n Qn an d therefore of physical 
significance: Whereas the original Hamiltonian of the Friedmann universe of scale factor 
exp Q driven by gravitons is the same as that of a Friedmann universe of the same scale 
factor driven by a massless minimally coupled scalar field, the transformed Hamiltonian is 
the same as that of a Friedmann universe with scale factor exp Q driven by a conformally 
coupled scalar or an electromagnetic field. A minimally coupled and a conformally coupled 
massless scalar field - or equivalently gravitational and electromagnetic radiation - driving 
a Friedmann background do not give rise to the same spacetime, beyond leading order in 
the Isaacson approximation. Only in the second case the effective stress-tensor is precisely 
traceless. To interpret a classical or quantum solution of the transformed system one must 
transform it back to variables expQ and q n (when the closed form of the solution is lost). 
The fact that the background spacetime with scale factor expQ is more than just one of 
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several equivalent choices of variable is obscured by the simplicity of our model: In a more 
complete model the physical scale factor could be determined with the help of any additional 
test matter field introducing a fundamental scale. The present paper derives solutions of 
the quantum system directly in the physical variables. This scheme is directly linked to the 
derivation of semiclassical gravity, while the solution of [EI] has the merit of being in closed 
form. 

The "Born-Oppenheimer" solution scheme might have been introduced by starting with 
the canonical transformation q n = (expQ)q n , Q = Q, giving rise to terms Ppq and pqpq 
in the transformed Hamiltonian. These are just the 2- and 4-particle interactions we have 
described. Seen this way, the 4-particle interactions should be neglected as being of quartic 
order in the perturbations. This means that we cannot consistently derive the equation of 
motion for o$ and higher order corrections in e. The author takes this as an indication that 
the physical content of a quantized perturbed minisuperspace model does not go beyond 
semiclassical gravity - higher order perturbative approximations to quantum gravity require 
a better approximation to the full classical action. 
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